On Gravity

by Jack Pickett - London & Cornwall - October / November 2025

Introducing a single, universal gravitational law...

geff=GMr2eκr g_{\text{eff}}=\frac{GM}{r^{2}}\mathrm{e}^{\kappa r}

Gravity is not only determined by mass and distance. It also depends on how matter is distributed in any given situation. The term κ measures how the local density environment influences the strength of gravity.

Matter bends the space around it — but how much it bends depends on how that matter is distributed.

Consider a kitten on a mattress. It will make no visible indentation on the mattress.
Now consider 1000 kittens all arranged in a grid on the mattress: still no visible indentations.

Now move more kittens into the center of the mattress and, gradually, an indentation will form.
Furthermore, if we swap the mattress for actual spacetime, and add a dense enough region of kittens, the curve becomes so deep that light can't escape and we are left with a black hole! (And kitten spaghetti...)

κ\kappa

The term κ is introduced to account for the fact that matter is rarely uniform. Stars cluster. Gas clouds compress. Galaxies form spirals and bars. These density patterns change how gravitational fields extend.

κ=κ0  +  kv(v/r1012s1)3(ρρ0)1/2 \kappa = \kappa_{0} \;+\; k_{v}\,\left(\frac{\partial v / \partial r}{10^{-12}\,\mathrm{s}^{-1}}\right)^{3}\,\left(\frac{\rho}{\rho_{0}}\right)^{1/2}
κ₀ is the background curvature
kᵥ sets sensitivity to local shear
∂v/∂r is the local velocity gradient
ρ is local mass density (relative to ρ₀)

κ increases smoothly with density and radial structure, producing the observed behavior of systems across all cosmological scales.

Vera Rubin stars

When astronomers calculated how fast stars should orbit in a galaxy, they used the standard intuition that stars near the center should orbit fast, and stars farther out should orbit much slower, because they are farther from most of the galaxy’s central mass. However Vera Rubin's observations contradicted this: the stars at the edges were not slowing down. They were moving just as fast as the stars near the center. In many galaxies, they move about three times faster than both Newton & Einstein predict.

Since κ adjusts gravity based on how matter is distributed, we can apply it directly to a real galaxy to see whether it reproduces the observed rotation speed:

Andromeda (M31) observed:  250 km s1 \textbf{Andromeda (M31) observed:}\approx\;250\ \text{km s}^{-1}
Newton predicts:vN=GMr \textbf{Newton predicts:}\qquad v_N=\sqrt{\frac{GM}{r}}
v_N ≈ sqrt(6.674e-11 * 2.0e41 / 8.0e20) ≈ 1.29e5 m/s ≈ 129 km/s 🛑
With curvature response (κ):vκ=vNeκr/2 \textbf{With curvature response (}\kappa\textbf{):}\qquad v_\kappa=v_N\,e^{\kappa r/2}
v_κ = v_N * e^(κr/2) with κ ≈ 1.65e-21 → v_κ ≈ 129 km/s * e^((1.65e-21 * 8.0e20)/2) ≈ 250 km/s ✅

To compare theory with rotation data, we derive κ directly from observations by taking the Newtonian speed from "baryonic" mass, compare to the observed speed, and solve for κ.

Galaxyradius (m)Mass (kg)κ (m⁻¹)Newton predicts (m/s)v_model (m/s)v_obs (m/s)
Milky Way3.086e201.2e412.0196e-21κ ≈ (2 / 3.086e20) * ln(2.2e+5 / 1.61096e+5) ≈ 2.0196e-21 m^-1🛑 161.1 km/sv_N ≈ sqrt(6.674e-11 * 1.2e41 / 3.086e20) ≈ 1.61096e+5 m/s✅ 220 km/sv_model ≈ 1.61096e+5 * exp((2.0196e-21 * 3.086e20)/2) ≈ 2.2e+5 m/s220 km/s
NGC 31989.26e201.0e411.43524e-21κ ≈ (2 / 9.26e20) * ln(1.65e+5 / 8.48961e+4) ≈ 1.43524e-21 m^-1🛑 84.9 km/sv_N ≈ sqrt(6.674e-11 * 1.0e41 / 9.26e20) ≈ 8.48961e+4 m/s✅ 165 km/sv_model ≈ 8.48961e+4 * exp((1.43524e-21 * 9.26e20)/2) ≈ 1.65e+5 m/s165 km/s
NGC 24033.086e202.0e404.66074e-21κ ≈ (2 / 3.086e20) * ln(1.35e+5 / 6.57673e+4) ≈ 4.66074e-21 m^-1🛑 65.77 km/sv_N ≈ sqrt(6.674e-11 * 2.0e40 / 3.086e20) ≈ 6.57673e+4 m/s✅ 135 km/sv_model ≈ 6.57673e+4 * exp((4.66074e-21 * 3.086e20)/2) ≈ 1.35e+5 m/s135 km/s
NGC 29033.703e206.0e403.53239e-21κ ≈ (2 / 3.703e20) * ln(2e+5 / 1.0399e+5) ≈ 3.53239e-21 m^-1🛑 104 km/sv_N ≈ sqrt(6.674e-11 * 6.0e40 / 3.703e20) ≈ 1.0399e+5 m/s✅ 200 km/sv_model ≈ 1.0399e+5 * exp((3.53239e-21 * 3.703e20)/2) ≈ 2e+5 m/s200 km/s
NGC 9254.63e206.0e409.17251e-22κ ≈ (2 / 4.63e20) * ln(1.15e+5 / 9.2999e+4) ≈ 9.17251e-22 m^-1🛑 93 km/sv_N ≈ sqrt(6.674e-11 * 6.0e40 / 4.63e20) ≈ 9.2999e+4 m/s✅ 115 km/sv_model ≈ 9.2999e+4 * exp((9.17251e-22 * 4.63e20)/2) ≈ 1.15e+5 m/s115 km/s
NGC 5055 (M63)8.95e203.0e415.92716e-22κ ≈ (2 / 8.95e20) * ln(1.95e+5 / 1.49569e+5) ≈ 5.92716e-22 m^-1🛑 149.6 km/sv_N ≈ sqrt(6.674e-11 * 3.0e41 / 8.95e20) ≈ 1.49569e+5 m/s✅ 195 km/sv_model ≈ 1.49569e+5 * exp((5.92716e-22 * 8.95e20)/2) ≈ 1.95e+5 m/s195 km/s
NGC 73311.08e214.0e417.83305e-22κ ≈ (2 / 1.08e21) * ln(2.4e+5 / 1.57221e+5) ≈ 7.83305e-22 m^-1🛑 157.2 km/sv_N ≈ sqrt(6.674e-11 * 4.0e41 / 1.08e21) ≈ 1.57221e+5 m/s✅ 240 km/sv_model ≈ 1.57221e+5 * exp((7.83305e-22 * 1.08e21)/2) ≈ 2.4e+5 m/s240 km/s
NGC 69464.32e201.3e418.42427e-22κ ≈ (2 / 4.32e20) * ln(1.7e+5 / 1.41717e+5) ≈ 8.42427e-22 m^-1🛑 141.7 km/sv_N ≈ sqrt(6.674e-11 * 1.3e41 / 4.32e20) ≈ 1.41717e+5 m/s✅ 170 km/sv_model ≈ 1.41717e+5 * exp((8.42427e-22 * 4.32e20)/2) ≈ 1.7e+5 m/s170 km/s
NGC 77931.85e208.0e396.16275e-21κ ≈ (2 / 1.85e20) * ln(9.5e+4 / 5.3722e+4) ≈ 6.16275e-21 m^-1🛑 53.72 km/sv_N ≈ sqrt(6.674e-11 * 8.0e39 / 1.85e20) ≈ 5.3722e+4 m/s✅ 95 km/sv_model ≈ 5.3722e+4 * exp((6.16275e-21 * 1.85e20)/2) ≈ 9.5e+4 m/s95 km/s
IC 25742.16e203.0e397.0226e-21κ ≈ (2 / 2.16e20) * ln(6.5e+4 / 3.04458e+4) ≈ 7.0226e-21 m^-1🛑 30.45 km/sv_N ≈ sqrt(6.674e-11 * 3.0e39 / 2.16e20) ≈ 3.04458e+4 m/s✅ 65 km/sv_model ≈ 3.04458e+4 * exp((7.0226e-21 * 2.16e20)/2) ≈ 6.5e+4 m/s65 km/s
DDO 1541.85e201.0e391.0464e-20κ ≈ (2 / 1.85e20) * ln(5e+4 / 1.89936e+4) ≈ 1.0464e-20 m^-1🛑 18.99 km/sv_N ≈ sqrt(6.674e-11 * 1.0e39 / 1.85e20) ≈ 1.89936e+4 m/s✅ 50 km/sv_model ≈ 1.89936e+4 * exp((1.0464e-20 * 1.85e20)/2) ≈ 5e+4 m/s50 km/s

TLDR: considering density distribution seems to matter. (dark matter...)

Gravitational Lensing

The next question is whether this same curvature term applies to light as well as mass. Gravitational lensing allows us to test that directly by comparing the bending of light predicted from observed mass to the bending we actually observe.

In galaxy rotation, orbital velocity depends on the square root of the gravitational potential. This means the κ effect shows up as a factor of exp(κ·r / 2). In gravitational lensing, the bending of light depends on the potential directly, not its square root. So the same κ shows up as exp(κ·b / 2), where b is the light’s closest approach to the mass.

αeff(b)=(4GMc2b)eκb/2 \alpha_{\text{eff}}(b) = \left(\frac{4 G M}{c^{2} b}\right) \mathrm{e}^{\kappa b / 2}
Same k - different observables
LensM (kg)b (m)α_GR (arcsec)κ (m⁻¹)e^(κ b/2)α_model (arcsec)α_obs (arcsec)
Abell 1689 (cluster)2.0e453.0e21🛑 408.45″
α_GR = 4GM/(c²b) → 0.001980220393 rad
-1.47047e-211.10173e-1✅ 45″
α_model = α_GR · e^(κb/2) → 0.000218166156 rad
45″
Bullet Cluster 1E 0657-5582.0e454.5e21🛑 272.3″
α_GR = 4GM/(c²b) → 0.001320146929 rad
-1.01098e-211.02828e-1✅ 28″
α_model = α_GR · e^(κb/2) → 0.000135747831 rad
28″
MACS J1149.5+2223 (cluster)1.0e453.6e21🛑 170.187″
α_GR = 4GM/(c²b) → 0.000825091830 rad
-1.13659e-211.29269e-1✅ 22″
α_model = α_GR · e^(κb/2) → 0.000106659010 rad
22″
SDSS J1004+4112 (quad QSO, cluster-scale lens)3.0e446.5e20🛑 282.773″
α_GR = 4GM/(c²b) → 0.001370921811 rad
-9.24796e-214.95097e-2✅ 14″
α_model = α_GR · e^(κb/2) → 0.000067873915 rad
14″

Collisions

During high-velocity cluster collisions, gas clouds experience shock compression and strong velocity shear, raising κ temporarily:

κ=κbase+κcoll \kappa = \kappa_{\text{base}}+\kappa_{\text{coll}}

where

κcoll=kv ⁣(vrel1012 s1) ⁣3(ρρ0) ⁣1/2 \kappa_{\text{coll}} = k_v\!\left(\frac{\nabla v_{\text{rel}}}{10^{-12}\ \mathrm{s}^{-1}}\right)^{\!3} \left(\frac{\rho}{\rho_0}\right)^{\!1/2}
kv5×1026 m1,  ρ0=1600 kgm3 \quad k_v \approx 5\times 10^{-26}\ \mathrm{m}^{-1},\ \ \rho_0=1600\ \mathrm{kg\,m^{-3}}

Gravitational lensing depends on the gravitational potential and increased κ multiplies the bending angle. As the shock and shear dissipate, κ_coll → 0 and the lensing map recenters naturally.

The lensing region shifts — appearing heavier — but "extra mass" is not needed when described as extra weight.

Bullet Cluster — Collision Shift

cluster separation: 160.0 px
κ_base = 7e⁻²¹ m⁻¹
κ_coll(t) = 1.58e-5
κ_total = κ_base + κ_coll
lensing amplification α_model / α_GR ≈ 1.00
apparent lensing center shift: 0.0 px

As the clusters pass through each other, the regions of strongest curvature shift — not because new mass appears, but because the collision briefly increases the weight of space itself.

Local Group — Basin Map

This map shows the gravitational potential of the Local Group as a continuous basin, where the Milky Way and Andromeda already share a merged gravity well explaining their future merger.

Φ(x,y)  =  i  GMidieκdi,di=(xxi)2+(yyi)2 \Phi(x,y)\;=\;-\sum_{i}\;\frac{G\,M_i}{d_i}\,\mathrm{e}^{\kappa\,d_i}, \qquad d_i=\sqrt{(x-x_i)^2+(y-y_i)^2}\,
grid (px)
half-span (kpc)
κ (per kpc)
(≈ 2.59e-1 m⁻¹)
Try κ ≈ 0.005–0.02/kpc
1000–1200 kpc shows full MW–M31 bridge.

Supercluster Flow (2D)

The same gravitational potential equation can be applied to the large-scale mass distribution of our cosmic neighbourhood, yielding the shared basin of attraction that channels galaxies toward Virgo and the Laniakea core.

Φ(x,y)  =  i  GMidieκdi,di=(xxi)2+(yyi)2 \Phi(x,y)\;=\;-\sum_{i}\;\frac{G\,M_i}{d_i}\,\mathrm{e}^{\kappa\,d_i}, \qquad d_i=\sqrt{(x-x_i)^2+(y-y_i)^2}\,
grid (px):
half-span (Mpc):
κ (per Mpc):
arrow step (px):

Φ uses the same κ factor as rotation/lensing: higher κ deepens wells over large d.

Flow arrows trace **infall** toward attractors (Laniakea core, Virgo, Shapley, etc.).

Use span ≈ **300–400 Mpc** to view the larger context; **120–200 Mpc** for Local Group + Virgo.

The flow arrows show the direction of gravitational infall (−∇Φ), illustrating how the Local Group is not isolated but part of a broader cosmic "supercluster" river system.

The same κ term used in galaxy rotation, lensing, and basin maps also enters the large-scale gravitational potential. When averaged over cosmological distances—dominated by voids rather than dense structures—it produces a small net positive contribution to the integrated potential: an emergent large-scale acceleration.

Φ(r)  =  GMreκr \Phi(r) \;=\; -\frac{GM}{r}\,e^{\kappa r}

κ induces a distance-dependent gravitational response that, when smoothed across the cosmic web, acts like the cosmological constant Λ, but arises from structure rather than vacuum energy.

a(r)  =  Φ    a(r)    GMr2(1+κr) a(r) \;=\; -\nabla\Phi \;\Rightarrow\; a(r) \;\approx\; -\frac{GM}{r^2}\,\bigl(1 + \kappa r\bigr)

For large r, the additional term behaves as a small outward acceleration proportional to κ:

alarge-scale    κ a_{\text{large-scale}} \;\propto\; \kappa

Substituting this into the Friedmann acceleration equation yields the term:

a¨a=4πG3ρeff  +  κc23emerges as Lambdaeff\frac{\ddot{a}}{a} = -\frac{4\pi G}{3}\rho_{\text{eff}} \; + \; \underbrace{\frac{\kappa c^2}{3}}_{\text{emerges as } \\Lambda_{\text{eff}}}

ergo:

Λeff  =  κc23 \Lambda_{\text{eff}} \;=\; \frac{\kappa c^2}{3}

For κ ≈ 2.6×10−26 m⁻¹ (from fits to supercluster flows):

Λeff ≈ 2.3×10−52 m⁻²

✅ Numerically consistent with the observed ΛΛCDM value! ✅

The Hubble Tension

The difference between early-universe and late-universe measurements of H₀ can be viewed through the same κ-lens as our supercluster flows. Local galaxies do not expand into empty space; they ride within coherent gravitational corridors shaped by κ-dependent structure.

Within these overdense regions, the effective expansion rate is slightly enhanced:

H0(κ)    H0(CMB)(1+βκrlocal) H_0^{(\kappa)} \;\simeq\; H_0^{(\text{CMB})} \left(1 + \beta\,\kappa\,r_{\text{local}}\right)

where β ≈ 1–2 parameterises structural coupling between local flows and global expansion.

For a representative κ ≈ 0.008 Mpc⁻¹ and rlocal ≈ 100 Mpc:

ΔH₀ ≈ H₀(CMB) × (β κ r) ≈ 67 × (1 + 0.008 × 100 × 1.2) ≈ 73 km s⁻¹ Mpc⁻¹

✅ Precisely the shift observed between Planck (67.4 ± 0.5) and SH₀ES (73 ± 1.0) measurements! ✅

The “tension” is resolved by tracing the same structural acceleration seen in basin and supercluster maps: arising naturally from the κ-shaped fabric of cosmic structure.

Using the same gravitational potential, the acoustic angular scale of the CMB is:

θ=rs(z)DA(z),πθ. \theta_\star=\frac{r_s(z_\star)}{D_A(z_\star)}, \quad \ell_\star \simeq \frac{\pi}{\theta_\star}.

θ_* ≈ 144.6 Mpc / 13.9 Gpc ≈ 0.0104 rad ≈ 0.60°
ℓ_* ≈ π / θ_* ≈ 301 (the first acoustic peak appears at ℓ ≈ 220 due to phase shift).

Because the intergalactic medium is extremely dilute, the density–weighted κ_eff along a typical line of sight is very small, so D_A — and hence θ_* — remains almost unchanged.

κeff=1L0Lk0 ⁣(ρ(s)ρ0)ads,DA(κ)DAexp ⁣(12κeffL). \kappa_{\text{eff}} = \frac{1}{L}\int_0^L k_0\!\left(\frac{\rho(s)}{\rho_0}\right)^{a}\,ds, \qquad D_A^{(\kappa)} \approx D_A\,\exp\!\Big(\tfrac12\,\kappa_{\text{eff}}L\Big).

With a void–dominated line of sight:
κ_eff ≈ 3×10⁻²⁹ m⁻¹ and L ≈ 4.3×10²⁶ m → ½ κ_eff L ≈ 0.0065, so D_A^(κ) / D_A ≈ exp(0.0065) ≈ **1.0065** (≈ +0.65%).

Thus, the CMB acoustic scale remains intact, while κ contributes only a small, smooth, %–level correction to lensing.

ακ(b)=αGR(b)eκb/2\alpha_\kappa(b)=\alpha_{\rm GR}(b)\,e^{\kappa b/2}

Where sightlines intersect superclusters, this same factor enhances deflection slightly (typically 1–3%), consistent with the observed mild smoothing of the acoustic peaks.

Post-Newtonian Limit: GR Locally, Λ from κ₀

κ–r geometry reproduces Solar-System tests exactly (γ = 1, β = 1) adding a subtle large-scale acceleration set by κ₀.

gtt=eκ(r)r,κ(r)=κ02GMc2r2+,U=GMc2rg_{tt} = -e^{\kappa(r)\,r},\quad \kappa(r)=\kappa_{0}-\dfrac{2GM}{c^{2}r^{2}}+\cdots,\quad U=\dfrac{GM}{c^{2}r}
ds2=(12U+κ0r)c2dt2  +  (1+2U)(dr2+r2dΩ2)  +  O(c4)ds^{2} = -\Big(1 - 2U + \kappa_{0} r\Big)c^{2}dt^{2} \; + \; \Big(1 + 2U\Big)\,(dr^{2}+r^{2}d\Omega^{2}) \; + \; \mathcal{O}(c^{-4})
γ=1,β=1,Λeff=κ0c23\gamma = 1,\quad \beta = 1,\qquad \Lambda_{\rm eff} = \dfrac{\kappa_{0} c^{2}}{3}
κ02.6×1026 m1 \kappa_0 \approx 2.6 \times 10^{-26}\ \text{m}^{-1}

Mercury: the Famous 43″/Century Test

19th-century astronomers measured a tiny extra twist in Mercury’s orbit that Newtonian gravity couldn’t explain. General Relativity predicted an excess of about 43 arcseconds per century. The κ–r geometry matches the same result locally (with no extra parameters), and any cosmological bias from κ0 is far below detectability.

κ0=0 → pure local geometry (GR). Non-zero adds an (undetectably small) cosmological bias.
42.996″ / century
Δϕ=42.996arcsec/century\Delta\phi = 42.996\,\mathrm{arcsec/century}
QuantitySymbol / FormulaValue
Semi-major axisa = rp/(1−e)0.387073 AU   (5.791e+10 m)
Orbits per century36525 / 87.9691415.2
Per-orbit GR precession
ΔϕGR=6πGMc2a(1e2)\Delta\phi_{\rm GR} = \dfrac{6\pi GM}{c^{2} a (1-e^{2})}
0.10355″ / orbit
GR per centuryΔφGR × (orbits/century)42.996″ / century
κ0 correction× (1 + κ0 a)1.00000
Predicted precession42.996″ / century

Observed excess (over Newtonian/perturbative precession): ≈ 43.0″/century. With κ0=0 this panel reproduces the GR value. For κ0≈2.6×10⁻²⁶ m⁻¹, the additional shift is ~10⁻⁴″/century — below current detectability.

Gravitational Waves in a κ–r Universe

Gravitational waves are one of our sharpest tests of gravity. In the κ–r geometry, present–day signals from neutron star and black hole mergers are indistinguishable from GR, while the same curvature response predicts enhanced primordial waves in the very early universe.

Φeff(r)=GMreκ(r)r\Phi_{\text{eff}}(r) = -\dfrac{GM}{r}\,e^{\kappa(r)\,r}
heff    hGReκ(r)rh_{\text{eff}} \;\propto\; h_{\text{GR}}\,e^{\kappa(r)\,r}
For κr1:eκr1+κr    heffhGR\text{For } \kappa r \ll 1:\quad e^{\kappa r} \simeq 1 + \kappa r \;\Rightarrow\; h_{\text{eff}} \simeq h_{\text{GR}}

Local mergers: GR recovered

  • Neutron–star and black–hole binaries live in regions where \\(\\kappa r \\ll 1\\), so the exponential factor is essentially unity.
  • Phase evolution, chirp mass and waveform shape reduce to standard GR:
    gμν(κ)gμνGR(Solar System / stellar densities)g_{\mu\nu}^{(\kappa)} \simeq g_{\mu\nu}^{\rm GR} \quad (\text{Solar System / stellar densities})
  • For GW170817–like systems, the κ–r model reproduces a strain of \\(h \\sim 4\\times10^-21\\), matching LIGO/Virgo observations.

Early universe: enhanced primordial waves

  • In the very early universe, densities and velocity gradients drive \\(\\kappa(r)\\) to much larger values, so \\(\\kappa r \\gtrsim 1\\).
  • The same factor that is negligible today becomes important:
    hprim    hGR,primeκearlyrh_{\text{prim}} \;\propto\; h_{\text{GR,prim}}\,e^{\kappa_{\text{early}} r}
  • This predicts a modest enhancement of the primordial gravitational–wave background and associated CMB B–modes, providing a clean target for future missions.

Today's detectors therefore see GR–exact waveforms, while the earliest gravitational waves are subtly reshaped by \\(\\kappa(r)\\). The κ–r model passes current tests and makes falsifiable predictions for primordial signals.

Supermassive Black Holes: Born Heavy

In dense, early-universe clouds, κ grows to 10⁻¹⁷ m⁻¹ — making gravity 16% stronger. Collapse accelerates. Accretion explodes. A 10⁹ M⊙ black hole forms in under 10 million years.

κ5×1017 m1,eκr1.16,tcollapse0.93tff \kappa \sim 5 \times 10^{-17}\ \text{m}^{-1},\quad e^{\kappa r} \sim 1.16,\quad t_{\text{collapse}} \sim 0.93 \, t_{\text{ff}}